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Abstract. We study the production of cosmic rays (CRs) in supernova remnants (SNRs), including the reaccelera- 
tion of background galactic cosmic rays (GCRs) — thus refining the early considerations by Blandford & Ostriker 
(1980) and Wandel et al. (1987) — and the effects of the nuclear spallation inside the sources (the SNRs). This 
combines for the first time nuclear spallation inside OR sources and in the diffuse interstellar medium, as well 
as reacceleration, with the injection and subsequent acceleration of suprathermal particles from the postshock 
thermal pool. Selfconsistent CR spectra are calculated on the basis of the nonlinear kinetic model. It is shown that 
GCR reacceleration and CR spallation produce a measurable effect at high energies, especially in the secondary 
to primary (s/p) ratio, making its energy-dependence substantially fiatter than predicted by the standard model. 
Quantitatively, the effect depends strongly upon the density of the surrounding circumstellar matter. GCR reac- 
celeration dominates secondary CR production for a low circumstellar density. It increases the expected s/p ratio 
substantially and flattens its spectrum to an almost energy-independent form for energies larger than 100 GeV/n 
if the supernovae explode on average into a hot dilute medium with hydrogen number density A^'h ~ 0.003 cm"''. 
The contribution of CR spallation inside SNRs to the s/p ratio increases with increasing circumstellar density 
and becomes dominant for A^h^I cm~^, leading at high energies to a flat s/p ratio which is only by a factor of 
three lower than in the case of the hot medium. Measurements of the boron to carbon ratio at energies above 
100 GeV/n could be used in comparison with the values predicted here as a consistency test for the supernova 
origin of the GCRs. 

Key words, theory - cosmic rays - shock acceleration - supernova remnants - reacceleration - secondary to 
primary ratio 



1. Introduction 
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It is a widely accepted hypothesis that supernova rem- 



Released from SNRs, CRs experience energy- 
dependent diffusion and, possibly, large scale convective 
transport in the interstellar medium (ISM). During 
their propagation in the Galaxy, the primary relativistic 
nuclei interact with interstellar gas nuclei and produce 
lighter secondary relativistic nuclei as a result of nuclear 
spallation. The secondary nuclei have approximately the 
same energy per nucleon as the parent primaries. 




the Galaxy. Supernovae occur randomly in the Galaxy (al- 
though with some correlation in space) and have a compli- 
cated time evolution. The selfconsistent nonlinear kinetic 
theory of diffusive shock acceleration nevertheless gives 
rather definite predictions on the shape, the amplitude 
and the time evolution of the spectra of CRs accelerated 
at the shock during the different phases of SNR evolution 
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Boron nuclei represent an example of pure secondaries 
and the Boron-to-Carbon (B/C) ratio is an example of the 
secondary-to-primary (s/p) ratio in CRs. The abundance 
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of secondary nuclei, measured over a wide energy range, 
yields important information on the process of CR trans- 
port in the Galaxy (e.g. Berezinskii et al. 11990(1. Up to 
approximately 100 GeV/n the measured power law spec- 
trum of secondary nuclei is steeper than the spectrum of 
primary nuclei. This leads to the conclusion that primary 
GCRs are mainly produced in compact sources (SNRs for 
example) and that their spectrum n(ek) oc ^"^ ' ^'^^ kinetic 
energies ek ^ 1 GeV/n, is steeper than the spectrum, 
A''(ek) oc e^'^" of the particles produced in the sources, 
the so-called source cosmic rays (SCRs): the GCR power 
law index 7g = 7s -I- /i is considerably larger than the 
SCR power law index 7s with a value /i r anging between 
At = 0.3 and 0.7. 

The diffusion of relativistic particles in the Galaxy, af- 
ter they have left the compact sources (assumed here to be 
SNRs) , may be accompanied by an additional distributed 
reacceleration in the turbulent ISM. This reacceleration 
can not be strong for particles with energies 1 GeV/n 
<ek^lOO GeV/n as evident from the observed energy- 
dependent decrease of the abundance of secondaries in this 
energy range (Hayakawa 196?). Some weak reacceleration 
is possible and may in fact explain the observed energy 
dependence of the s/p ratios with a characteristic peak 
at about 1 GeV/n (see Jones et al. I2U01I and references 
therein). In the "minimal" model, the stochastic reaccel- 
eration is produced by the same randomly moving inho- 
mogeneities which are assumed to be responsible for the 
spatial diffusion of CRs in the galactic magnetic field. The 
velocity of these inhomogeneities is close to the Alfven ve- 
locity Ca ~ 30 km/s. The characteristic time of distributed 
reacceleration is approximately equal to the escape time 
from the Galaxy for particles with magnetic rigidities of 
about 1 GV (the magnetic rigidity is i? = pc/Ze , where p 
is the particle momentum, Ze denotes the particle charge, 
and c is the speed of light). Escape dominates over reac- 
celeration at high rigidities i? >> 1 GV, so that the effect 
of distributed reacceleration in the interstellar turbulence 
is weak at high energies. In detail the importance of dis- 
tributed reacceleration for the formation of the GCR spec- 
trum remains unclear. It may well be that the peaks in 
the s/p ratios are caused by a specific dependence of the 
CR diffusion on energy or by the effect of large scale con- 
vection on the transport of low energy CRs (Jones et al. 
I2()01|l . It should be stressed however that the distributed 
reacceleration through interstellar Alfvenic turbulence is 
negligible at energies ek^20 GeV/n in all versions of the 
galactic CR propagation. 

The scenario described is valid for the major part 
of the GCRs which left the source regions and do no 
more interact with SNR shocks propagating in the ISM. 
However, a wandering energetic particle has a finite proba- 
bility to again meet a SNR shock which efficiently acceler- 
ates CRs, and therefore to undergo strong reacceleration. 
The background primary and secondary CR nuclei experi- 
ence strong reacceleration and significant transformation 
of their energy spectra after interaction with high velocity 



shocks. Strongly reaccelerated primaries and secondaries 
at the shock have flat spectra typical of the source spec- 
trum. These particles are mixed with and diluted by a 
large amount of background GCRs after their release from 
their parent SNR. In principle they may become recogniz- 
able above some energy as a flat component of galactic 
secondary nuclei. 

The effect described above was studied by Wandel et 
al. (1987) (see also the original paper by Blandford & 
Ostriker 1980|l who found that a flattening of the s/p ratio 
can be expected at high energies due to reacceleration by 
strong SNR shocks. At the same time, applied to an actual 
galactic ISM, this model leads to the conclusion (Wandel 
1988) that one has to expect a reacceleration effect mainly 
due to old and weak SNR shocks of Mach numbers M < 2 
which also essentially influence the CR spectra at GeV 
energies and make the CR mean escape length from the 
Galaxy as small as 3 to 5 g cm~^. These conclusions were 
based on a simplified plane wave approach. All the SN 
shocks were considered to be of the same strength and the 
final predictions depend substantially on this suggested 
strength. 

We note also that it is hard to believe that such weak 
shocks are even able to accelerate relativistic particles in 
the same manner as a plane shock. Indeed, to be efficiently 
accelerated by a spherical shock of size Rg and speed Vg, 
particles should have a sufficiently small diffusion coeffi- 
cient (e.g. Berezhko I199()f) 

K<0.1RsVs. (1) 

Shocks of size Rs ~ 30 pc and speed Vg « 300 km/s 
(WandelEnHHI) give the restriction k<3 x 10^^ cm^/s. Such 
a small value can only be produced near the shock front 
by selfexcited Alfven waves, because the CR diffusion co- 
efficient due to interstellar turbulence exceeds this value 
even at GeV energies (e.g. Berezinskii et al. ll990|l . On the 
other hand the spectrum of particles accelerated by such 
weak shocks is so steep, n oc p~'^ with 7 « 4, that, even 
if extended into the relativistic energy range p > mc, it 
contains too small an amount of energy to substantially 
amplify the ISM turbulence. This means that efhcient ac- 
celeration of relativistic particles in SNRs stops earlier, at 
higher Mach numbers Af « 4, when the shock compres- 
sion ratio drops significantly below the asymptotic non- 
relativistic value cr = 4, and freshly accelerated particles 
become unable to provide a high level of turbulence near 
the shock front. Therefore it is assumed that at this age 
t — TsN, SNRs release all previously accelerated CRs into 
the surrounding ISM. 

The energy increase of relativistic CRs upon encoun- 
tering the SNR of age t > Tsn is negligible because their 
relatively large diffusive length / = k/Vs > 0.1i?s reduces 
their energy gain at the shock front almost to the adi- 
abatic energy loss in the expanding downstream region. 
Therefore we neglect the acceleration and/or reaccelera- 
tion of CRs by these old SNRs. 

Since according to the above physical arguments the 
maximum age and the corresponding number of SNRs 
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which effectively accelerate CRs is presumably much 
smaller than assumed by Wandel H1988|l . the effect of reac- 
celeration of GCRs is also much smaller. Nevertheless it 
is still important, especially at high energies for the case 
of a low-density ISM. 

In this paper we shall in addition take the effect of nu- 
clear spallation inside the sources into account. The energy 
spectrum of these source secondaries is harder than that 
of reaccelerated secondaries. Therefore it plays dominant 
role at high energies for a high-density ISM. 

Note that the CR source spectrum N{e]t} which is the 
resultant overall spectrum of CRs accelerated during the 
whole active period of SNR evolution depends on the GCR 
spectrum n(ek) since not only the injection of suprather- 
mal particles leads to the formation of N{e]^) but the 
reacceleration of existing background GCRs contributes 
as well. These spectra iV(ek) should therefore be deter- 
mined through a selfconsistent solution of the equations 
which describe the GCR propagation in the Galactic con- 
finement volume with all accompanying effects, together 
with the equations which describe the SNR evolution and 
the CR acceleration. Fortunately this rather complicated 
problem can be simplified without loss of accuracy. Indeed, 
due to the steep GCR spectrum at relativistic energies 
ek ^ 1 GeV/n only the injection of GCRs with energies 
Ck^l GeV/n plays a role for the formation of Af(ek)- Since 
the GCR spectra n(ek) at these subrelativistic energies are 
known from experiment, we use them as an physical input 
quantities, like other ISM parameters, in our selfconsistent 
nonlinear kinetic model for CR acceleration in SNR, so 
that -/V(ek) can be considered as consistently determined 
together with the GCR spectra n(ek)- 

We also want to emphasize that the shape of the overall 
CR spectrum N{ek,t) differs at each stage of SNR evolu- 
tion significantly from the spectrum of freshly accelerated 
CRs n{ek,t), especially during the late Sedov phase. In 
this late phase, when the shock is relatively weak and thus 
CR backreaction not significant, the freshly produced CR 
spectrum has a pure power law form n oc p~^, bounded 
by an exponential cutoff at some maximum momentum 
Pm{t)- At the same time, the overall CR spectrum is in- 
fluenced by the nonlinear shock modification during the 
previous active period of the SNR evolution, by the adi- 
abatic cooling of CRs in the downstream region, and by 
the so-called escape effect (see Berezhko et al. '1996' for 
details) and therefore has a much more complicated form. 
Only at very low energies ^ 1 GeV it has a form simi- 
lar to n{p,t), whereas at higher energies it becomes pro- 
gressively flatter and is bounded by the cutoff momentum 
Pmax = Pmito) > Pm{t), a valuc which is determined by 
the very early Scdov phase fBerezhko I1996|l . This is an 
additional and significant difference between our model, 
based on the kinetic nonlinear time-dependent approach, 
and Wandel's ri988) estimates, based on a plane- wave 
approximation, which does not include these important 
physical effects. 

In section 2 we present the relevant relations within 
the standard leaky box model. In section 3 the GCR in- 



jection problem within the nonlinear kinetic model for 
CR production in SNRs is described and the relevance 
of secondary nuclei production in SNRs for the resultant 
source spectra -/V(ek) and for the s/p ratio are estimated. 
Numerical results for the boron and carbon spectra pro- 
duced in SNR and the expected B/C ratio are presented 
and discussed in section 4. The final section contains the 
main conclusions. 

2. Source contribution to the secondary to 
primary ratio 

We shall use the simple leaky box model to describe the 
transport and nuclear fragmentation of comic rays in the 
Galaxy. It is well known that the more adequate diffusion 
model of cosmic ray propagation with an extended flat 
halo gives essentially the same results as the leaky box 
model for the abundances of the stable primary and sec- 
ondary nuclei. The approximate equivalence of these two 
models holds provided that cosmic ray nuclear fragmen- 
tation and reacceleration occur in regions much thinner 
than the cosmic ray halo size. This last condition is ful- 
filled in the Galaxy where the disk thickness of the order 
0.5 kpc is much smaller than the halo thickness (which is 
at least an order of magnitude larger) (Berezinskii et al. 
11990). 

For a simple estimate we consider the case when the 
secondary nuclei of a species "s" are the result of spalla- 
tion of some primary CR nuclei of species " p" , and the 
contribution of all other primary parent nuclei can be ne- 
glected. Within this approach the differential number den- 
sity of some primary CR species 7ip(ek) is determined by 
the balance equation 

n N 

— = -rj-V's.N ~ (TpNgVUp, (2) 
Tc Vc 

where Vc is the residence volume occupied by the GCRs, 
Tc is the mean residence time in this volume Vc, cTp is 
the spallation cross section, 7Vp(ek) is the total differential 
number of CRs created during the entire evolution of a 
single SNR (overall CR spectrum), i^sn is the galactic SN 
explosion rate, v is the particle speed, A'g — pa/nip is 
the ISM number density, po is the ISM density, and mp 
denotes the proton mass. 

Within the same approach the number density of sec- 
ondary nuclei Hs obeys the equation 

— = tf'^sn + cTpsNgVUp - a^NgVUs, (3) 

Tc Vc 

where Ups is the partial cross section for creation of sec- 
ondary nuclei in hadronic collisions of primaries with the 
ISM gas nuclei, Ng is the total number of secondary CRs 
created during the evolution of a single SNR. Compared 
with the standard leaky box model the above equation 
contains an additional source term which describes the 
production of secondaries in SNRs. In the above equa- 
tions all variables are functions of the same kinetic energy 
per nucleus Ck- 
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One can easily find the s/p ratio from the above two 
equations: 



1 + a.NgVTc Np (1 + asNgVT^) 



(4) 



Since the residence time Tc in the GCR confinement vol- 
ume is a decreasing function of CR energy for ek ^ 
1 GeV/n, we have 



aNo-VTc < 1, 



(5) 



and the expression for the ratio can be drastically simpli- 
fied: 



n 



where we have introduced the escape length 



X = mpNgVTc, 



(6) 



(7) 



in the form of the mean matter thickness traversed by 
GCRs in the course of their random walk in the Galaxy. 
Expression jnj shows that apart from the usual term pro- 
portional to the escape length x, the s/p ratio contains an 
additional term which describes the source contribution 
in the production of secondaries. Since the primary spec- 
tra Np{e\^), and the secondary spectra 7Vs(ek), presumably 
produced in SNRs, have a similar energy dependence, we 
expect the term Ng/Np to dominate at sufficiently high 
energies in the s/p ratio. 

One can also get an expression for the escape length x 
from Eq. in terms of the observed s/p ratio rg = Us/up 
and the source s/p ratio rg = Ns/Np-. 



mpirg-rg) 



(8) 



In the simplest case when the source produces only 
primary CRs (rg = 0), we have the usual expression 



X = mp/(crps/rg - CTs), 



(9) 



which makes it possible to extract the value of x from 
the measured s/p ratio rg. Note that in the general case, 
when an unknown number of secondary CRs are produced 
in the CR sources, there is no one to one correspondence 
between rg and x. 

3. Acceleration and reacceleration of CRs in SNRs 

There are two different suprathernial particle populations 
in the ISM which are injected into the diffusive shock ac- 
celeration process in SNRs. The first and most general one 
is the injection of some fraction of the postshock thermal 
particle distribution. It occurs for all ions present in the 
background medium and usually supplies enough particles 
to convert an significant part of the SN shock energy into 
that of an energetic particle population. 

The second possibility is the acceleration of pre- 
existing GCR particles which have a sufficiently high en- 
ergy Ek^lOO MeV/n so that they participate naturally 



into the acceleration process. Note that for those elements 
which are strongly underabundant in the thermal ISM, 
like Li, Be, B, this is the only practical possibility for 
acceleration. To distinguish these two different injection 
mechanism we use here the term "acceleration" for the 
first case and "reacceleration" for the second. 

We shall employ here a simple CR injection model, in 
which a small fraction rj of the incoming thermal protons 
is instantly injected at the gas subshock with a speed that 
exceeds the postshock gas sound speed Cs2 by a factor 
A > 1 (Berezhko et al lH^ : 



(10) 



Here iVg = p/mp is the gas number density, and the sub- 
scripts 1(2) refer to the point just ahead (behind) the 
shock. 

GCR nuclei (primary and secondary) have a spectrum 
which in the subrelativistic region rises as a function of 
Ek due to ionization and nuclear losses. The spectrum has 
a maximum at ek = cgcr ~ 600 MeV/n and then falls 
of in a steep power law for higher energies. Therefore in 
the case of reacceleration it is assumed that the existing 
GCR population is injected at the SN shock front into the 
diffusive acceleration with this energy 



Pinj = PGCR, Mi 



GCR- 



(11) 



Here ^Vqcr is the total number of GCR nuclei per unit 
volume and pgcr is their mean momentum, that corre- 
sponds to eccR- 

Since the primary CRs come from both injection mech- 
anisms, we can estimate their relative role for the final 
CR production in SNRs. In order to do that we com- 
pare the overall CR spectra (i.e. the particle numbers dN 
in the momentum interval dp) N{p) and N'^'^ij)) which 
correspond to these two injection mechanisms. The over- 
all spectrum of CRs reaccelerated in each single SNR at 
P > Pgcr is determined by the simple relation 



VsnA^gcr / P 



Pgcr Vpgcr 



(12) 



where Vsn is the volume of the SNR at the latest evo- 
lutionary stage where efficient CR acceleration can occur. 
For the sake of simplicity we assume that the overall spec- 
trum of CRs produced during the active period of SNR 
evolution has a pure power law form with index 7s. 

After their release from the parent SNRs these reac- 
celerated CRs occupy the confinement volume Vc more or 
less homogeneously with a number density Uq^^ that can 
be found from the leaky box equation (O (we neglect here 
the collision term): 



'GCR 



/Tc = iVVsN/K. 



(13) 



Since GCRs have a power law spectrum tigcr p 
their total number density is 



NgCR = ?^GCR(PGCR)^'GCR/(7g " 1) 



(14) 
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and therefore we can write 



"GCR 
'^GCR 



(7g - l)Vc Vpgcr 



7g-7= 



(15) 



Note that this relation is valid for primary as well 
as for secondary elements. The only difference between 
them is the different momentum dependence of the ratio 
'^GCr/'^GCR- For primary elements as it follows from the 
equation Q 



7g = 7s + 



(16) 



where power law index fj, determines the momentum de- 
pendence of the residence time Tc oc p^^. Therefore the 
ratio J^GCPi /^gcr is independent of p in this case. 

The spectrum of secondary elements, produced in the 
Galactic disk according to the equation Q , is considerably 
steeper, because 



7g = 7s + 2m, 



(17) 



and therefore in this case the contribution of GCR reac- 
celeration 



"-gcr/"gcr oc (p/pgcr)'' 



(18) 



increases with energy. 

The overall CR spectrum is formed during the active 
period of SNR evolution which lasts up to the time when 
the SN shock becomes too weak to accelerate efficiently 
a new portion of freshly injected particles. According to 
calculations presented below this final stage corresponds 
to a shock size Rs ~ 10i?o in the case of an ISM with 
hydrogen number density A^h = 1 cm~'^, and Rs ~ 3i?o 
in the case of a hot ISM with — 0.003 cm^'^, where 
Rq = (3A/oj/47rpo)^/'^ is so the called sweep-up radius, 
Afoj is the ejecta mass, and po — N^mp = 1.4iVHTOp is 
the ISM density. Taking into account that Tc/Vc = Tg/Vg, 
where Tg is the residence time in the galactic disk vol- 
ume Vg = 2.5 X lO'^*' cm^, and i/sn = 1/30 yr~\ we 
have for energies about 5 GeV, where Tg — 4.6 x 10^ yr, 
a value riQ^^/nQCB. = 0.04 for iVn = 1 cm^'^, and 
"^gcr/^-gcr — 0.3 for iVn = 0.003 cni^-^. One can see that 
in an ISM which has a number density TVh = 1 cm~'^ typ- 
ical for the galactic disk, the primary GCR reacceleration 
is not important. Only in the case of a rarefied ISM, like a 
hot ISM with a number density — 0.003 cm^'^, one can 
expect an important contribution of GCR reacceleration 
to their overall spectrum nGCR(ek)- The same conclusion 
has been reached for CR acceleration in the ISM of ellip- 
tical galaxies that has similar parameters (Dorfi & Volk 

Due to the additional factor {p/pgcr}^ the contribu- 
tion of reacceleration becomes dominant in the secondary 
GCR spectra at high energies. The value of the critical 
energy above which the secondary spectra are expected 
to be as hard as the primary ones depends upon the 
mean gas density TV^^^ inside the SNRs: according to 
the above formulae, it is expected to be 100 GeV/n for 



cm , and an order of magnitude larger if 
7V|CR _ I cm-3. 

There is an additional mechanism of secondary GCR 
production inside SNRs: primary nuclei like GCRs in the 
Galactic disk produce light secondary as a result of nu- 
clear spallation due to their interaction with the back- 
ground gas. One can easily estimate the relative signifi- 
cance of this mechanism compared with the previous one. 
The total production rate of secondaries in SNRs due to 
the spallation of primaries is 



-iSCR 



(19) 



where A^sn is the number of existing SNRs in the Galaxy 
and N^'^^ is mean number density of the thermal gas in 
SNRs. GCRs generate the same kind of secondaries with 
the rate 



GCR ^ N^CRy^np, 



(20) 



where N^'^^ 



1 cm^'^ is mean gas number density in the 
Galactic disk. Taking into account that the SNR number 
A^SN — JsNt'SN is determined by the CR confinement time 
TsN inside a SNR, and that the GCR number density rip 
is related to the overall CR spectrum produced in SNR 
A^p by the relation rip = A^pt^sNTg/V^, we have 



gSCR ^SCR 



SN 



IGCR 



A^gGCR Tg 



(21) 



Like reacceleration, this mechanism produces a substan- 
tially harder spectrum of secondaries compared with the 
spectrum created by the GCRs in the Galactic disk. It has 
the opposite dependence on the ISM density: it is more ef- 
ficient in a denser ISM. 

Since TgiM ~ 10^ yr and since Tg « 5 x 10® yr 
at GeV energies (e.g. Berezinskii et al. 11990(1 . we have 
(5^c:R/gGCR ^ 2 X 10"^ if on average SNe explode into 
the ISM with a gas number density N^'^^ = N^'^^ . This 
value is somewhat lower than the corresponding contri- 
bution of reacceleration to the production of secondary 
nuclei. At the same time, the spectrum of secondaries 
spectrum Ns{p) has in this case exactly the same form as 
that of the primaries, Np{p). As demonstrated below, in 
the case of reacceleration the ratio A^g(j>)/Afp(p) decreases 
slightly with momentum p due to the fact that secondaries 
are more effectively produced at the late SNR evolution- 
ary phase when the SN shock becomes weaker, whereas 
high energy primaries are mainly accelerated at the begin- 
ning of the Sedov phase, when the SN shock is extremely 
strong. Therefore one can expect that the secondary nu- 
clei production due to the hadronic SCR collisions with 
the gas nuclei makes a considerable contribution at high 
energies ek ^ 100 GeV/n. 

The production of secondaries due to spallation of pri- 
mary SCRs in SNRs is described in our model by the 
source term 



(?s(r,p, t) = apsvNgfp{r,p, t) 



(22) 
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in the transport equation for the distribution function 
fs{r,p,t) of secondary nuclei. Here fp{r,p,t) is the distri- 
bution function of primary nuclei and Ng is the local gas 
number density, both of which are also functions of SNR 
evolutionary time and position. The distribution function 
f{r,p, t), which determines all relevant CR characteristics, 
is calculated in our model selfconsistently for all CR ele- 
ments considered (see Berezhko et al. 119961 for details). 

The actual situation with the CR spectra in SNRs is 
rather more complicated than suggested by the above es- 
timates. Due to nonlinear acceleration effects, CR spec- 
tra produced in SNRs are not pure power laws. Besides 
that particles with different energies are effectively pro- 
duced during different phases of SNR evolution, whereas 
the above simple estimates suggest in fact that the en- 
tire overall CR spectrum N(jj) is produced at an epoch 
characterized by the SNR volume Vsn and age Tsn- To 
reach more definitive conclusions about the role of CR 
reacceleration in SNRs, we need to perform a selfconsis- 
tent consideration. This can only be done numerically. 

Note that the source CR spectra -/V(ek) which are the 
time integral spectra of the CRs accelerated during the 
whole active period of SNR evolution, depend on the GCR 
spectra ri(ek) since not only the injection of suprathcrmal 
particles leads to the formation of -/V(ek) but the injection 
and the reacceleration of existing background GCRs con- 
tribute as well. This means that the spectra N{£]^) should 
be determined as the result of a selfconsistent solution of 
the equations which describe the GCR propagation and 
nuclear transformation in the Galactic confinement vol- 
ume, together with the equations, which govern SNR evo- 
lution and CR acceleration. Fortunately this rather com- 
plicated problem can be considerably simplified without 
loss of accuracy. Indeed, as discussed before, due to the 
steep GCR spectra at relativistic energies ek ^ 1 GeV/n, 
only the injection of GCRs with energies ek;$l GeV/n 
plays a role for the formation of the resultant spectra 
-/V(ek). Since the GCR spectra n(ek) at these subrelativis- 
tic energies are known from experiment, we use them as 
physical input like other ISM parameters in our selfconsis- 
tent nonlinear kinetic model for CR acceleration in SNR, 
so that the output of this model, which is the overall CR 
spectrum N{e\^), can be considered as consistently deter- 
mined together with the GCR spectra n(ek)- 

4. Results and discussion 

In order to study the effect of GCR reacceleration in SNRs 
on the s/p ratio we consider the elements boron (B) and 
carbon (C) because carbon is the most abundant element 
in the GCRs and plays a major role in the production of 
boron which is usually considered as a pure secondary. We 
have performed selfconsistent calculations of CR acceler- 
ation in SNRs, based on the kinetic nonlinear model, for 
the simple case of a uniform ISM with different densities. 

We use the values i?sN = 10^^ erg for the explosion 
energy and Mcj — 1.4Mq for the ejecta mass which are 
typical for SNe la in a uniform ISM. Note that the main 



o 


1 

7 / 


1 


T^, K 
2.2x10^ 


t/to : 

200 


ro 


7/ ■ X 




8.6x10^ 


500 ■ 








3.8x10^ 


1000 


to 










cn 
CN 

T 


//' 

- / 

7 


c/ 




■ 


in 

T 










00 

T 


1 


1 




\ \ " 

w 




-2 





2 


4 6 



log(p/mpC) 

Fig. 1. The overall momentum spectra of B and C nu- 
clei produced during the SNR evolution in the warm ISM 
(iVH = 0.3 cm-3, T = 10^ K, to = 367 yr) up to three 
different moments of time. The carbon spectra, starting 
from p/rripC ^ 1, are due to the acceleration of injected 
suprathcrmal particles and GCR reacceleration, whereas 
the boron spectra are the result of reacceleration of GCRs 
and the spallation of carbon and oxygen nuclei inside SNR. 
For clarity, the spectrum of purely reaccelerated GCR car- 
bon nuclei, corresponding to t/to = 1000, is shown by the 
thin solid line. 

fraction of the core collapse SNe has relatively small ini- 
tial progenitor star masses between 8 and 15 Mq which 
therefore do not significantly modify the surrounding ISM 
through the main sequence wind of the progenitor star 
(e.g. Abbott 1982 ). SNR evolution in this case is very sim- 
ilar to that of SNe la. 

Only in the case of type lb and type II SNe with mas- 
sive progenitor stars Mi > IBMq, SNRs expand into a 
nonuniform circumstellar medium strongly modified by 
the intense progenitor star wind. The SNR evolution looks 
differently compared to the case of a uniform ISM. At the 
same time, the main amount of CRs is expected to be 
produced when the SN shock propagates through the hot 
rarefied bubble, where the evolution and CR acceleration 
are roughly similar to that in the case of a uniform ISM 
which has the same parameters as the bubble (Berezhko 
and Volk 2000J. Therefore the contribution of SNRs ex- 
panding in the strongly modified circumstellar medium 
effectively increase the filling factor of the hot phase of 
ISM. 

We adopt an injection rate of suprathcrmal particles, 
characterized by the injection parameters = 10~* and 
A = 4, which are expected for a typical SNR (Volk et al. 
I2()()3(l . GCR nuclei B and C with number densities N-q = 
7.9 X IQ-^'^ cm-3 and Nc = 2.6 x IQ-^^ cm'^, respec- 
tively, are injected at a kinetic energy einj = 0.6 GeV/n 
which corresponds to the mean GCR energy for these el- 
ements. We consider three essentially different phases of 
the ISM: a diluted, hot ISM with hydrogen number den- 
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sity A^H = 0.003 cm^^ and temperature Tq = 10*^ K, a 
warm ISM with A^h = 0.3 cm^^ ^^j^^j j^^ ^ ^O" K, and an 
"average" ISM with = 1 cm^^ and To = 10"* K. The 
ISM magnetic field values Bp = 3 /iG, 5 /iG and 10 ^iG 
where taken for these three cases, respectively. 

Numerical results, which correspond to the so-called 
warm ISM with N-a = 0.3 cm~^, are presented in Figs.^ 
and|21 The overall momentum spectra of boron and carbon 

/•oo 

7V(p,t) = 167rV / drr^f{r,p,t), (23) 
Jo 

which include all particles accelerated during the SNR evo- 
lution up to time t, are shown in Fig.^for three successive 
times t. One can see that reacceleration of GCRs makes 
quite a small contribution to the overall spectrum of car- 
bon, in agreement with the earlier estimate. Note also that 
the late evolutionary phases are not very important for 
the highest energy part of the spectrum of primary nuclei: 
the C-spectrum at p > lO^WpC is formed mainly at the 
beginning of the Sedov phase {t ^ to), and it is almost 
unchanged at late phases t > lOO^o. There are two physi- 
cal factors which determine the above situation. The first 
is a geometrical factor which results in the most efficient 
acceleration of the highest energy CRs at the beginning 
of the Sedov phase fBerezhko .l996;l . Due to the decrease 
of the product RsVg, the maximum energy of freshly in- 
jected/accelerated CRs decreases with time in the Sedov 
phase. Therefore the late Sedov phases do not contribute 
significantly to the highest energy part of the resulting CR 
spectrum. 

The second factor is the decrease of the injected 
suprathermal particle momentum pinj oc Vg due to the 
shock deceleration. It diminishes the contribution of the 
late phases to CR production. As a result, the overall spec- 
trum of CR primary nuclei is mainly formed during the 
early Sedov phase when the shock is extremely strong, and 
therefore it is quite hard. 

In the case of secondaries the situation is very differ- 
ent. In the B-spectra two different components are clearly 
visible. The first one, peaking at p = pgcr ~ lOmpC, is 
due to the injection of GCRs with momentum pgcr ■ Since 
Pgcr does not depend on time, the efficiency of their accel- 
eration (reacceleration) progressively increases with time, 
roughly proportional to the total number of particles par- 
ticipating in the acceleration. This number is proportional 
to the shock volume Vsn- Since in this case the late SN 
evolutionary phases, when the SN shock is no more so 
strong as during the early Sedov phase, are important for 
secondary CR production, their overall spectrum is much 
steeper than the spectrum of primary nuclei (see Fig. ^ . 

The second component in the B-spectra is due to the 
spallation of carbon and oxygen nuclei. Its spectrum is 
substantially harder than that of the first component, be- 
cause it is directly related to the spectra of primaries. 
Therefore it contributes strongly at very high momenta 
p^lO'^nipC and also dominates at very low momenta 
P < PGCR: where B production due to reacceleration is 
negligible. 
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Fig. 2. The B/C ratio Vg — r^B/f^c as a function of ki- 
netic energy per nucleus. Solid, dashed and dash-dotted 
lines represent the ratio Vs = N-q/Nq for the case pre- 
sented in Fig. Dotted lines represent the fit of the ex- 
perimental data within the leaky box model (thin dots) 
on the one hand, and within the diffusive model with 
distributed reacceleration in the ISM {heavy dots) on the 
other. Experimental data (Engelmann et al. 11990(1 are also 
shown. 

On account of the above factors the ratio rg, which is 
exclusively due to the production of secondaries in SNRs, 
becomes steeper and increases with time, as one can see 
from Fig. [3 Therefore the effect of GCR reacceleration 
on the s/p ratio depends on the confinement time Tsn 
which corresponds to the age of the SNR when the shock 
becomes an inefficient accelerator and release of previously 
accelerated CRs into the Galactic volume proceeds. 

There are at least two conditions which determine the 
confinement time Tsn. 

The first one is the decrease of the shock Mach num- 
ber during the SNR evolution. When it becomes so low, at 
some stage t = t^, such that the shock compression ratio 
a drops below 4, the acceleration of freshly injected parti- 
cles becomes inefficient. This leads to the decrease of the 
turbulence level near the shock front and to subsequent 
escape of previously accelerated CRs due to increase of 
their diffusive mobility. This factor plays the major role 
in the case of a hot ISM. 

In fact, the so-called effective compression ratio 

def = a(l-l/A4), (24) 

rather than cr, plays the main role in CR acceleration and 
their final spectrum. Due to the fact, that Alfven waves, 
excited by the accelerated CRs upstream of the shock, 
propagate in upstream direction with Alfven speed Ca, the 
effective compression ratio of scattering centers seen by 
CRs CTef is always smaller than a. This effect is described 
in the above expression by the Alfven Mach number Afa — 
Cn/Vg. Therefore t^ is the time when cjcf drops below 4. 

The acceleration process may also stop at some stage 
t = te when the postshock temperature drops below 10^ K, 
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because radiative SNR cooling sets in strongly. Since the 
SNR looses a large amount of its energy, efficient CR ac- 
celeration terminates presumably at this stage. 

We adopt here a confinement time, which is the mini- 
mum of these two: 

TsN = min{t4, te}. (25) 

In the case of a warm ISM, we have a postshock gas tem- 
perature T2 = 1.3 X 10® K at t = 300^0, and therefore the 
confinement time is TgN = ^6 ~ 10'^ yr {to = 367 yr in this 
case). The relevant ratio rg in Fig. [21 is therefore given by 
the dashed line. 

The significance of GCR reacceleration can be under- 
stood if one compares the expected source s/p ratio rg 
with the ratio rg, extracted from the measured B/C ratio 
in the standard model. 

Since the boron nuclei are produced in the ISM not 
only by spallation of carbon but also by spallation of oxy- 
gen, the expression has to be rewritten in the form: 

UB X f "o\ , iVe (mp + (Tea;) 

nc trip + (Tbx\ nc J Nc(mp + aBX) 

This expression is used in order to determine the ex- 
pected B/C ratio with accounts for boron production in- 
side the GCR sources (in SNRs). In the case when the 
secondaries are not produced in CR sources we have from 
this expression 

^ KB+aoB)Xg^ 

nip + aBXg 

According to the standard leaky box model, the value 
of the escape length Xg « 14 g/cm^ at GeV energies with 
energy dependence Xg oc e^°'^ at ek > 5 GeV/n fits the 
existing B/C data (Engelmann et al. ll99Ufl . 

Note however that the value of the escape length de- 
rived from the same experimental data depends also on 
the GCR propagation model, especially at high energies 
where the existing experimental data do not constrain it. 

The galactic model with distributed stochastic reac- 
celeration predicts a different (weaker) decrease of the 
B/C ratio at high energies ek > 30 GeV/n than does 
the standard leaky box model. According to Jones et al. 
mm . Xg oc Ck model without reacceleration, 

and Xg oc fjj ^'^ in the model with stochastic reacceleration. 
Both models give the same s/p ratio at ek < 30 GeV/n. 

In order to demonstrate how significantly secondary 
CR production in SNRs changes the expected s/p ratio, 
we consider here two extreme predictions for Xg oc ej^^. 
The first one, with /i = 0.6, corresponds to the standard 
leaky box model (Engelmann et al. 1990). The second, 
with fj, = 0.3, corresponds to the diffusion model with 
distributed stochastic reacceleration (Jones et al. I2()()l(l . 
The B/C ratio calculated with these two functions a;g(ek) 
according to the expression (|26|l with Nb = 0, are shown 
in Figs.|2HSlby the dotted lines. The cross section values 
and the ratio no/«Cj required to calculate the B/C ratio, 
are taken from the paper by Engelmann et al. I)199()|l . 
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Fig. 3. The same as in Fig. |5|but for the case of a hot 
ISM (iVn ^ 0.003 cm-^, to = 1674 yr). 

In Figs. [21Blwe also present the ratio rg — Nb/Nq 
of the boron and carbon spectra, produced in SNRs by 
the combination of all processes considered up to dif- 
ferent times t/tQ and calculated for three different ISM 
densities. One can see from Fig. 12 which corresponds to 
A^H — 0.3 cm^'^, that the source ratio rg exceeds the B/C 
ratio Tg predicted by the leaky box model at an energy 
ek > 300 GeV/n and only at ek > 10 TeV/n it becomes 
larger than rg which corresponds to the diffusive model 
with distributed reacceleration if the CR confinement time 
is as large as Tsn = 4 x 10^ yr. Note however, that at 
such a late time t = lO^to = 4 x 10^ yr, the postshock 
temperature is already T2 « 4 x 10^ K, and in addition 
the effective compression ratio drops up to a^f = 3.25. 
Therefore the CR confinement time is taken to be equal 
TsN = 500^0 = 2 X 10^ yr, because the postshock tem- 
perature at this stage is T2 = 9 x 10^ K and the effective 
compression ratio has a value ad = 3.8 which is only 
slightly below 4. 

In Fig. 131 we present the same as in Fig. [21 but for the 
rarefied, so-called hot ISM which is characterized by a 
temperature Tq — 10® K and a hydrogen number density 
A^H = 0.003 cm^"^. One can see that the ratio Vg exceeds at 
energies ek^lOO GeV/n the measured value of B/C ratio 
already at t = 3 x 10"* yr. 

We note that already at this stage the effective com- 
pression ratio drops below 4 and the acceleration of 
freshly injected particles becomes progressively less effi- 
cient. Therefore one can assume that in this case Tsn = 
^4 « 3 x 10^ yr. Even for this relatively short confinement 
time the contribution of the reacceleration process in the 
B/C ratio is strong at energies ek greater than 10 GeV/n 
(solid curve in Fig.lSJ. Note that at such a low ISM den- 
sity the production of secondary nuclei in SNRs is entirely 
due to their reacceleration. 

At larger ISM densities secondary nuclei production 
due to the primary CR spallation becomes more impor- 
tant. This is illustrated in Fig.0I where the expected B/C 
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Fig. 4. The same as in Fig. |21but for the case of an ISM 
hydrogen number density A^h = 1 cm^"^ [t^ = 246 yr). 

ratio Ts is represented for A'h — 1 cm^^. It is seen that 
at high energies ek/^lC GeV/n the expected ratio Ts is 
higher than in the case iVn = 0.3 cm~'^. Therefore, we 
conclude that the lowest efficiency of secondary nuclei pro- 
duction in SNRs occurs in the case of ISM number density 
TVh = 0.3 cm~'^. It becomes larger both in a denser and a 
more diluted ISM because of efficient primary spallation 
and GCR reacceleration, respectively. Since at i = SOOfo 
the postshock temperature T2 = 9 x 10^ K is already lower 
than 10^ K, the CR confinement time is Tsn = ^6 ~ 10^ yr 
in this case even though the effective shock compression 
ratio (Tcf = 4.2 is still slightly higher than 4. 

Since the typical ISM density in which the majority 
of SN explosions occurs is not well known, we present in 
Figs. 13 and the possible range of the expected overall 
B/C ratios, calculated according to expression (|26() . The 
source ratios — N-q/Nq are taken from the above cal- 
culations and the CR escape length 

X = 5xg (28) 

was assumed to be some fraction (5 < 1 of the escape length 



determined from the relation (|27l) . The cross section values 
and the ratio no/nc, required to calculate the B/C ratio, 
are taken from the paper by Engelmann et al. II 199011 . To 
fit the existing data one needs 6 — 0.3, 0.7 and 0.9 for the 
ISM number densities Nn = 0.003 cm^^, 0.3 cm"'^ and 
1 cm""^, respectively. The B/C ratio is shown for these 
three densities. The production of secondary CRs in SNRs 
leads to a substantial decrease of the escape length x only 
if all SNRs explode into a diluted ISM with hydrogen num- 
ber density A^h — 0.003 cm^'^. In the two other cases, 
which are expected to be more probable, the relative de- 
crease of X is not so important. It is significantly lower 
compared with the estimates of Wandel 1)198811 . because 
we believe that the active period of SNR evolution stops 
much earlier than he assumed. 



0.1 1 10 100 1000 10'* 10^ 

GeV/n 

Fig. 5. The expected overall B/C ratio as a function of 
kinetic energy per nucleus. The three curves, correspond- 
ing to the three circumstellar hydrogen number densities 
A^H = 0.003 cm"3 /j^g)^ ^ q.S cm^^ (^dashed 

line), and iVu = 1 cm''^ {dash-dotted line) encompass a 
range given by the dotted area. 



It is seen that for A'h — 0.003 cm~'^ the boron produc- 
tion in SNRs leads to considerably higher values of the 
B/C ratio at energies ek > 100 GeV/n than the obser- 
vations show. There are however some experimental in- 
dications that the measured B/C ratio becomes indeed 
flatter at energies near 100 GeV/n. It is interesting to see 
that the theoretical curves do not change monotonically 
with gas density. In fact a weighted mixture of SN sites 
with very low densities, Nu = 0.003 cm-^ (hot ISM), and 
moderate densities A^h = 0.3 cm"'^ (warm ISM) is able to 
approximate the observations. More precise measurements 
at these energies are needed, in order to draw definite con- 
clusions about the actual role of GCR reacceleration. 

We note that the wide spread of the B/C ratio at high 
energies is due to the unknown correlations of SNR sites 
with gas density. Nevertheless it is clear from Figs. 
that at any given density the energy dependence of the 
B/C ratio is expected to be much flatter than predicted 
by the leaky box model. 

We have a similar situation in the case when the GCR 
escape length a;g(ek), predicted by the diffusive model with 
distributed reacceleration, is used for comparison with the 
theory, as one can see from Fig. IHl for the same three 
densities. Since the energy dependence of the ratio 7'g(ek) 
is much flatter in this case, the relative increase of the 
B/C ratio due secondary CR production in SNRs is much 
smaller than in the previous case. Nevertheless, the ex- 
pected value of the B/C ratio for ek<Jl TeV/n is still 
noticeably higher, by a factor of two on average, than pre- 
dicted by the model which neglects secondary CR produc- 
tion in SNRs. The CR escape length values are determined 
by the factor S = 0.2, 0.6 and 0.8 for A^h = 0.003 cm^^, 
0.3 cm"'^ and 1 cm^"^, respectively. 
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Fig. 6. The same as in Fig. |5l but comparing with the 
mean CR escape length Xg(ek) predicted by the diffusion 
model with distributed CR reacceleration. 

One could expect that we deal with a more realistic 
situation when a variety of existing SNRs in the Galaxy 
are expanding into ISM regions with different densities. 
In order to make a definite prediction about the expected 
B/C ratio one should determine properly weighted boron 
and carbon spectra: 

NsA^k) - ^AvlMNB,c{^k,Nl)/vsN (30) 

for a given SN explosion rate ^v^gj^i in an ISM with num- 
ber density 7V^. Unfortunately the explosion rates ^v^sn 
are not known. At the same time, as it is clear from Figs. 
5 and 6, the uncertainty in the B/C ratio is determined 
mainly by the value A^'g^, which corresponds to the hot 
ISM. Since the boron spectrum is much more sensitive to 
the ISM density than the carbon spectrum, we can ap- 
proximately write 

B/C = {B/C)h^v''sN/vsN+{B/C)^{l^^v''^N/vsN), (31) 

taking into account that within the ISM density range 
from 0.3 to 1 cm~'^ the ratio B/C changes not so sig- 
nificantly. Here {B/C)h and {B/C)w are the B/C ratios 
produced in hot and warm ISM, respectively. This formula 
makes it possible to estimate the expected B/C ratio for 
a given value Ai/gj^ /i/sn without making additional cal- 
culations. 

5. Summary 

Our considerations based on the selfconsistent kinetic non- 
linear model for CR acceleration in SNRs show that the 
reacceleration of the existing GCRs and the spallation of 
primary CRs due to their collisions with the gas nuclei in 
SNRs strongly influence the energy spectra of secondary 
elements, like Li, Be, B. Due to this additional mechanism 
the spectra of secondaries become significantly flatter at 
high energies ek > 100 GeV/n. This effect can be directly 



studied from the s/p ratio: at high energies where reaccel- 
eration is important, the s/p ratio should be flatter than 
at lower energies l;$ek;S100 GeV/n. 

Compared with previous considerations based on the 
plane- wave approach fWandel [l988|l . our model selfconsis- 
tently includes a number of physical factors: i) nonlinear 
shock modiflcation due to the CR backreaction; ii) CR adi- 
abatic cooling inside the expanding SNR interior; iii) the 
so-called CR escape effect. In addition, all stages of SNR 
evolution contribute to the resulting overall CR spectra, 
in contrast to the assumption by Wandel that only very 
late SNR phases play a role. 

Qualitatively, the efficiency of secondary CR nuclei 
production in SNRs depends in an important way on the 
density of the ISM in which SNRs are exploding. The lower 
the ISM density is, the larger is the volume of the SNR 
which is reached during the active period of its evolution. 
Correspondingly larger numbers of background GCRs are 
swept-up and reaccelerated compared with the number of 
particles injected from the thermal postshock distribution. 
If SNRs are on average exploding into a diluted ISM, then 
the influence of reacceleration becomes significant already 
at Ek ^ 10 GcV/n. This is in contradiction to the existing 
data for the B/C ratio and means that the typical ISM 
density "seen" by SNe is greater than 0.003 cm^'^. 

An increase of the ISM density leads to a decrease of 
the reacceleration effect. At the same time the production 
of secondaries in SNRs due to the spallation of primaries 
becomes more important. If the typical density is between 
about A^H = 0.3 and 1 cm"-^, reacceleration and spallation 
produce roughly equal contributions to the B/C ratio at 
Ek ~ 1 TeV/n, whereas at higher energies spallation be- 
comes dominant. Quantitatively the effect depends upon 
the mean GCR escape length and its energy dependence. 
In the case Ah = 0.3 cm^"^ boron production in SNRs be- 
comes dominant for ek > 0.5 TeV/n within the leaky box 
model, and comparable within the diffusion model with 
distributed reacceleration for ek > 3 TeV/n. 

The value of the CR conflnement time Tsn inside 
SNRs plays the crucial role for the secondary CR pro- 
duction inside SNRs: the larger Tsn is, the greater is the 
number of secondaries that are produced. In the case of 
a diluted ISM with number density Nh — 0.003 cm^"^, 
CR conflnement is presumably restricted to times smaller 
than when the shock becomes too weak to support a 
high level of MHD turbulence near the shock front due 
to the CR streaming instability. In the case of a much 
denser ISM with Nh ^ 1 cm"'^ the active period of CR 
production in SNRs is assumed to terminate when the 
postshock temperature drops below 10® K. For intermedi- 
ate ISM densities Nh ^ 0.3 cm~^ both effects are equally 
important. Due to these two physical factors the CR con- 
flnement time Tsn is significantly lower in our model com- 
pared with previous considerations (Wandel 1988), where 
these factors were ignored. This is one of the reasons why 
the effect of secondary CR production in SNRs is signif- 
icantly lower in our model compared with the prediction 
of Wandel [MHI 
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If on average SNe explode in the ISM with A^h = 
1 cm^'^ boron production in SNRs leads to an ex- 
tremely flat, almost energy- independent B/C ratio at e^, > 
1 TeV/n and at > 5 TeV/n in the two above cases, re- 
spectively. A detection of this flat energy-dependent ratio 
would be consistent with GCR production in SNRs. 

The actual mean galactic escape length x is smaller 
than in the models which neglect the production sec- 
ondary nuclei inside the sources and only take into ac- 
count spallation during the random walk of primaries in 
the Galactic disk after release from their sources. The ex- 
tent of this reduction depends upon the mean ISM den- 
sity at the SNR sites: x is from 10% to 30% smaller 
if SNRs on average explode in an ISM with hydrogen 
number density TVjj between 0.3 and 1 cm~'^, respec- 
tively. This effect is much smaller compared with the es- 
timates of Wandel H1988|l , who suggested efficient CR ac- 
celeration/reacceleration by much older and weaker SNR 
shocks. Only if the ISM number density is as low as 
A'h = 0.003 cm^'^ one expects the mean escape length 
to be smaller by a factor of 3 compared with the standard 
models. 
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